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Abstract. We present a model independent approach to FFs, which al- 
lows to determine uniquely the non-singlet combination {D^—Dd)^^^'^ 
of the kaon FFs, independently in the three different processes of kaon 
inclusive production: e^e" annihilation, pp-coUisions and IN semi- 
inclusive DIS. Only C and SU(2) invariances are used and no assump- 
tions on FFs and PDFs, the result holds in any order in QCD. The non- 
singlets obtained in this way, would allow to test the commonly used 
assumptions in global fit analysis and to test universality. This method 
is applied to kaon production data in + -I- X, = K^,K^ 

and {Du — Dd)^*'^^ is obtained without any assumptions. This re- 
sult is compared to the existing parametrizations, obtained in global 
fit analysis. 



1 Introduction 

In order to extract correct information about the quark and lepton interactions from 
high energy experiments, the precise knowledge of two basic quantities are required - 
the parton distribution functions (PDFs), that describe the distribution of the nucleon 
momentum among its constituents - the quarks and gluons, and the fragmentation 
functions (FFs) , that describe the distribution of the parton momentum among the 
hadrons, which the parton fragments into. Both the PDFs and the FFs are non- 
perturbative objects, to be inferred from experiment. Quantum Chromodynamics 
(QCD) does not provide a definite picture for their calculation, it determines only the 
order by order calculation in perturbation theory for their energy scale dependence, 
described by the Q^-evolution DGLAP equations. 

While PDFs are relatively well determined, there are lot of uncertainties about the 
FFs at present. The importance of the FFs has become evident only the last decade, 
with the start of the new generation of high-energy experiments with a final hadron h 
detected like semi-inclusive deep inelastic scattering (SIDIS) and hadron production 
in pp-collisions. The above information becomes still more viable for the LHC, where 
detecting different final hadrons will be our window for the New Physics expected. 
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At present different parametrizations about the FFs, obtained from global fit 
to data, exist |1I2I3I4] . Three characteristic features about them hold: 1) they use 
different theoretical assumptions about unfavoured FFs, 2) they describe the data 
and 3) they don't agree among themselves. 

Recently, the problem with the kaon FFs has become most appealing when the 
COMPASS collaboration [S] presented two very different numbers for the strange- 
quark polarization As = J dx As{x, Q^) using the same data but two different parametriza- 
tions for the FFs: DSS [1^ and EMC :4;: 

As = -0.01 (DSS), As = -0.04 (EMC). (1) 



2 The difference cross sections 

2.1 Cross sections with K"^ and 

We consider charged and neutral kaon production in the three semi-inclusive pro- 
cesses: 

e+ + e- ^ K + X (2) 
l + N ^l + K + X (3) 
p + p^ K + X (4) 

where K stands for either charged or neutral kaons, K — , , . Note that 
measuring neutral kaons does not introduce new FFs in these cross sections as SU(2) 
invariance of strong interactions relates neutral and charged kaon FFs: 

T~)K++K- _ j~)K°+K'> T^K+ + K- _ T^K° + K° t^K++K- _ j~,K°+K° 

We show that the measuring the difference of charged and neutral kaons: 



in any of the three processes (HI) - (jH), always measures the same NS combination 

{Du — Dd)^^^^ ■ This result is obtained without any assumptions about the PDFs 
and FFs, using only the SU(2) relations (O for the kaons and, being a consequence 
of a symmetry, holds in any QCD order. The results for the three processes are [B]: 

1) For the z-distribution in e+e^ — > (7, Z) ^ K + X we have: 

daf:^t^' -2A1^ ^ g ( „2 _ ,2)(i _^ g ^ ) ^K^A- (^^ q2) 

where ctq — 47rag^/3 s and Bq are the electroweak charges of the quark g, Cq are the 
perturbatively calculable Wilson coefficients, z is the fraction of the momentum of 
the fragmenting parton transferred to the hadron h: z — /E, where E^ and E are 
the CM energies of the final hadron and the initial lepton, and ^/s = 2E. 

2) The considered difference cross-sections ([3]) for SIDIS on proton and deutron 
targets are given by: 



C^^--'"^'" = l[iu + d) « (1 + I^Cqq) 2g .g Cqq] « D^^.^^- 
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Here q = u + u and Cab are the perturbatively calculable Wilson coefficients. 
3) For pp^ K + X, K = K^, K° we obtain: 



X (da^^f + da^J ~ dai^ ~ rfaff ) . (10) 

Here the sum over a, b is over all partons that contribute, fa^b sue the PDFs, and 
da-^-^ {s, t, u) are the partonic cross sections for the inclusive processes a + b ^ c + X 
that contribute. They are functions of the corresponding Mandelstam variables: 

S={Pa+Pb)'^, t={pa-PcT, U={pb-PcT 
Pa = XaPA, Pb = XbPs, Pc = /z. (11) 

(In NLO s, t and u are independent variables.) In addition, that the expression for 
the difference cross section tr^ is considerably simpler than the corresponding cross 
sections for single kaon production - in NLO it is only 8 partonic cross sections that 
contribute to cr^, versus 21 to cr^ . 

Eqs. ([7|- p^ are four independent measurements that determine uniquely the 
NS combination of the kaon FF D^^+^~ . As D^*+^~ is a NS, no new FFs will 
enter through its evolution. As these expressions use only SU(2) invariance, the 
verification of any of them might serve as a test for SU(2) invariance of the kaon FFs, 
used in all analysis. 

In addition, eqs. ([7])- © allow to compare the NS obtained in e+e" at rather high 
~ m^, with those from SIDIS at quite low . This would provide a challenging 
test of the Q^-evolution and universality of the FFs. 

As these expressions are model independent, it would be interesting to compare 
the resulting NS to the NS obtained with the existing parametrizations extracted 
from e^e~ data, obtained with various assumptions. This we shall do in Section [3l 

2.2 Cross sections with 

Let's consider charged kaon production in: 

l + N ^l + K + X (12) 
p + p^K + X, (13) 

K = K^,K~. We show that, using only C-invariance, the difference cross sections 
a'^*-^' = a^'^ - <7^' in (HD) and ^ determine uniquely the NSs D^*-^~ and 
~^ without any assumptions. The expressions in NLO are 0: 



a. 



a 



~ [Auv D^^-''' +dv® Df-""'] + a,Cg,), (14) 

K 

d 

K 

PP 



^ ~ [{uv + dv) ® (4i?„ + i?d)^^-^n ® (1 + asCqq) (15) 



a^p-'^ c^[Lu®uv®Du+Ld®dv®Dd\'^ . (16) 
Here Lq are known functions of the unpolarized PDFs and partonic cross sections: 
Lu{x, t, u) = u{x) dS{s, t, u) + [d{x) + S{x)] da'^^f, (s, t, u) + g{x) da'^^l^^^^ {s, t, u) 
Ld{x, t, u) = d{x) dt:{s, t, u) + [u{x) + s{x)] dal^, (s, i, u) + g{x) da^^-'i'^^ (s, i, u) 
Ls{x, t, u) — s{x) dX!{s, t, u) + [u{x) + d{x)] da'^^, {s, t, u) + g{x) dd^^g'^^'^ {s, t, u) 
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where dS is the combination 
dlJ 



da^^f{s,t,u) + ^dal'i-1^''{s,t,u) 



q = q + q. 



(17) 



If we assume D, 



— 0, used in all global analysis, from ([Til - ([TS]) we obtain: 



a. 



p 

K+-K- 
-'d 

PP 



{uv + dv) <E) (l + as Cqq) ® D!^^-^~ 
dz 

dxadxb— [Lu{xb,t,u) (g) Uv{Xa) + Lu{Xa,t, u) Uv{xb)] ® D. 



(18) 
(19) 



(20) 



Eqs. (dHJ - (PO)) present 3 independent measurements for the NS D^^^^ , obtained 

with the only assumption = and all a^^~^ are fitted with just one FF. 

As ~^ — has been the only assumption, each of these equations would be a 
test of just this assumption. 



3 /^f_7^" from e+e data 

On basis of eq. Q, analyzing the available data on and if^-production in e~^e~- 
annihilation, we extracted for the first time D^_J'^ without any correlations to 
other FFs. There are several arguments for choosing e+e~ data: 1) these are the 
most accurate and precise data, 2) the NS component doesn't contain unresummed 
soft gluons divergences, which allows to use lower values of 2: (z > 0.001) than in 
global fit analysis (z > 0.1, DSS [I]; z > 0.05, HKNS [S] and AKK08 0), 3) for 
the NS component a next-next to leading order (NNLO) fit is possible as both the 
splitting and coefficient functions are known to NNLO, while in global fit analysis 
only NLO calculations are available at present. 

The data used in our analysis are in the energy range -^/s = 12 — 189 GeV, but 
not all data appear equally important. Eq.© implies that the sensitivity of cr'^ to 

^u~d^ is determined by the s-dependence of (e^ — e^)(s). This difference is the 
biggest away from the Z-pole, i.e. most important for our studies are the data with 
small s, y/s < 60 GeV and big s, ^/s > 110 GeV. Unfortunately, the data around the 
.^-pole, where the most abundant and precise data exist, almost does not contribute. 

Our analysis would be easy if we had data on and K'^ production at the same 
z and s. However, as this is not the case, we proceed in three steps [7]: 1) we combine 
the measurements into 7 energy intervals: y/s = 12-14.8; 21.5-22; 29-35; 42.6-44; 58; 
91.2; 183-189 GeV, 2) in each interval we fit data purely phenomenologically and 

3) we parametrize D^_^^ at a factorization scale /io — \/2 GeV: 

D'^l+^-{z,^,l) = nz%l - z)" + n'z"'(l - z)y (21) 

and fit simultaneously the parameters in the parametrizations for and D^ J~^ 
in a combined analysis of the data on and K^. In our analysis we keep the hadron 
mass m/j as a fitting parameter. 

Small z- values imply small hadron energies, which may become of the order or 
smaller than the mass of the hadron. That's why including the small-z and/or small-s 
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data the effects of the hadron mass, i.e. rrih ^ 0, become significant and have to be 
taken into account. In our analysis we must distinguish kinematically between the 
experimentally measured fractions of the final hadron energy E^, and the modulus of 
its momentum jp^^j to ^/s: 

These quantities are equal and equal to the light cone scaling variable 
only at m/j = 0, x = z in LO in QCD. 

The NS is presented on Fig. [B left. On the same figure D^^+^' , 

obtained from the global fits of DSS, AKK08 and HKNS, are presented as well. The 
differences between our result for D^_^^ and those of global fit analysis might be 
due mainly to 1) the different assumptions - in DSS and HKNS Df* = D^* = 
D^^ = D^^ is assumed; in AKK08 = D^* is used; we make no assumptions; 

2) the small z-data, for the first time included. In order to get insight on the effect 

of the low-z data, on Fig [U right we present D^^J'^ with different cuts on the 
lower values of z. The result shows that data at z > 0.01 is not enough to impose 
constraints on the NS. 

At the end we point to a striking feature of our result - the negative value for 

^u-d^ at about z < 0.3. In order to check this result and find out its origin we 
performed our analysis with different parametrizations for D^_^^ and this does 
not seem a genuine feature of the NS - choosing D^^^^ = nz°'{l — z)*" in stead of 
(PT|) yields ffj^/^" > 0, D^^+^~ = appears also compatible with the data. The 
quality of the fit only slightly worsens: from Xdof = = 2-2 for D^_/^ < 0, to 

Xdof — 2.4 for positive and zero values of the NS, N — 730 is the number of points 
in the fit. However, in all cases the quality of the fit is worse than the one obtained 
in the purely phenomenological fit to - Xdof = l-l- 

Also the value of the fitted kaon mass, obtained in the purely phenomenological fit 
to K'^ - rriK = 320 MeV, differs significantly from the one obtained in the combined 

and analysis based on the QCD formula ([7]) - rriK = 124 MeV. (The true kaon 
mass is niK — 494 MeV.) A possible reason for this big difference might be that other 
low-Y^ and small- a; effects, such as higher twists and mass effects of resonances from 
which kaons may be produced, that have not been accounted for in our analysis, can 
also be absorbed in the fitted value for vhk- The performed analysis are very sensitive 
to hadron mass effects and these effects strongly affect our fits. More accurate data 
will be needed to determine how important these other effects are. 

Our analysis shows that the non-singlet D^_^ can be determined uniquely in a 
combined fit to e~^e~ data, but more precise data is needed to fix it better. 
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Fig. 1. up: D^^^'^ in this paper at NLO with resumed logarithms [label AC], and from 
DSS, HKNS and AKK08 sets; down: the NS in this paper with different cuts on the data 
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